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Abstract 

We investigate the sparticle spectrum in models of gauge-mediated supersymmetry 
breaking. In these models, supersymmetry is spontaneously broken at an energy scale only 
a few orders of magnitude above the electroweak scale. The breakdown of supersymmetry 
is communicated to the standard model particles and their superpartners by "messenger" 
fields through their ordinary gauge interactions. We study the effects of a messenger sector 
in which the supersymmetry-violating _F-term contributions to messenger scalar masses are 
comparable to the supersymmetry-preserving ones. We also argue that it is not particularly 
natural to restrict attention to models in which the messenger fields lie in complete SU (5) 
GUT multiplets, and we identify a much larger class of viable models. Remarkably, however, 
we find that the superpartner mass parameters in these models are still subject to many 
significant contraints. 
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1. Introduction 



The masses of the superpartners of the Standard Model (SM) particles should not greatly 
exceed the TeV scale if supersymmetry is to solve the hierarchy problem associated with the 
ratio Mz/Mpianck- However, this fact by itself tells us surprisingly little about the scale AsusY 
at which supersymmetry is ultimately broken. It is also necessary to have an understanding of 
the mechanism by which supersymmetry breaking is communicated from its original source to 
the fields of the Minimal Supersymmetric Standard Model (MSSM). If gravitational or other 
Planck-suppressed interactions communicate supersymmetry breaking, then AsusY is perhaps 
10^^ GeV or so. While this scenario has received the most attention in the last decade, it 
is hardly inevitable. Another possibility [||, |2| is that the ordinary gauge interactions are 
responsible for communicating supersymmetry breaking to the MSSM through their couplings 
to a messenger sector of chiral superfields, which in turn couple directly or indirectly to the 
fields which break supersymmetry. 

In the "minimal" model of gauge mediated supersymmetry breaking (GMSB) j^, all of 
the soft supersymmetry-breaking interactions of the MSSM are determined by just a few free 
parameters. Perhaps the most attractive feature of this type of model is that the masses 
generated for squarks and sleptons with the same SU{3)c x SU{2)l x U{1)y quantum numbers 
are automatically degenerate, so that flavor-changing neutral currents are suppressed without 
additional assumptions. This feature depends only on the fact that ordinary gauge interactions 
are flavor-blind, and will be true in a much larger class of models than just the minimal GMSB 
model. 

This class of models has another feature which may allow it to be dramatically confirmed at 
existing or currently planned collider facilities. Because local supersymmetry is spontaneously 
broken at a relatively low scale, the lightest supersymmetric particle is the gravitino (the spin 
3/2 superpartner of the graviton), with a mass that is entirely irrelevant for collider kinematics 
(but not for cosmology |^]). The next-to-lightest supersymmetric particle (NLSP) can therefore 
decay into its SM partner and the gravitino. In the case that the lightest neutralino (A'^i) is 
the NLSP, one has the interesting decay|Q, ^ iVi ^ 76 as long as the photino content of A^i 
is non-zero. The decay length for this process depends on the ultimate scale of supersymmetry 
breaking AsusY, according to 



r(iVi ^ 7G) 




(1) 
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where = \NiiCOs9w + A''i2 sin 6'vi/ P (in the notation of is the photino content of Ni. 
Since in GMSB the typical i^'-term responsible for supersymmetry breaking can correspond 
to AsusY of order 10^ or 10^ TeV, it is quite possible that this decay can occur (at least a 
significant fraction of the time) inside a typical detector, with many interesting phenomenolog- 
ical consequences [Q, ^, |^ ^, ^, 0, |ll], 12, 13, 14]. If it is sufficiently heavy A''! can also have 



decays into ZG and hG, with decay widths which suffer, however, from very strong kinematic 



suppression [12 1 . 



Recently it was pointed out ^, ^ that a single 6677 + event observed at CDF could 
be naturally explained^ by GMSB (and other theories with a light gravitino). This event had 
an energetic electron and positron, two energetic photons each with pseudorapidity \rj\ < 1 and 
transverse energy Et > 30 GeV, and large missing transverse energy > 50 GeV. The SM 
and detector backgrounds for such events are reputed to be extremely small. This event can be 
explained by GMSB as either selectron pair production or chargino pair production, but only if 
iVi is the NLSP, and if AsusY is less than about 10^ TeV. These are not automatic consequences 
of all models, and therefore will give (if taken seriously, which clearly should not be considered 
mandatory!) non-trivial theoretical constraints. 

Moreover, the discovery signatures of supersymmetry with a prompt decay A'^i — > 7G are 

so spectacular that it is possible to set quite strong bounds even with existing Tevatron data. 

In contrast to the usual supersymmetry search strategies, one can obtain a very high detection 

efficiency at the Tevatron for the inclusive signal 77 + X + with suitable cuts on the 

transverse energy and isolation of the photons, and on the total missing transverse energy. In 

|12| it was argued that with the present 100 pb^^ of data at the Tevatron, it should be possible 

to exclude a lightest chargino (Ci) mass up to 125 GeV and neutralino masses up to about 70 

GeV, assuming gaugino mass "unification" relations as in the minimal GMSB model. In this 

paper we will discuss other models which do not share this feature. Even when all assumptions 

about gaugino mass relations are abandoned, however, it was argued in that one can still 

find a model independent bound m^_^ > 100 GeV as long as > 50 GeV (to supply energetic 

photons) by exploiting the inclusive 77 + X + signal. These bounds are quite competitive 

with and somewhat complementary to what can be done at LEP upgrades. However, it should 

be kept in mind that these bounds all assume that the decay A'^i 'jG occurs within the 

detector 100% of the time. This is not necessary, even to explain the CDF 6677 + event, 

^The event can also be explained in the usual MSSM framework without a light gravitino, if parameters are 
chosen so that the radiative 1-loop decay N2 Nij dominates^, The parameter space in which this can 
occur will be largely but not entirely explored at LEP161 and LEP190. 
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which only requires that some non-negUgible fraction of A'^i decays occur within the detector. 
If most decays occur outside the detector, then one would expect many more single photon 
events than diphoton events, with unfortunately a much larger SM background, and much 
more difficult challenges for simulation studies. Thus for example the discovery mode at LEP2 
from e~^e~ NiNi could be predominantly jf! rather than 77-^. We should also note that in 
a significant fraction of the models to be studied in this paper, A^i cannot be the NLSP anyway 
unless it is higgsino-like. 

While the minimal model of GMSB is quite elegant and can explain the CDF 6677 + 
event, it is important to consider what all the related alternatives might be, especially in 
setting discovery and exclusion strategies. Future phenomenological studies should therefore 
take into account the full richness of model-building possibilities, which undoubtedly extend 
far beyond the minimal GMSB model and in sever al different directions 0, ||, 1^, ||, |2|. 



In this paper we will begin to explore a few such possibilities. In section 2 we develop the 
formalism for arbitrary messenger sector field content including the effects of arbitrary masses 
(from scalar VEVs and i^-term breaking) in the messenger sector. In section 3 we will examine 
the discrete model space allowed by generalizing the particle content of the messenger sector 
to include possibilities which do not form complete GUT multiplets. We will argue that it 
is not particularly unnatural or even inelegant to consider such generalizations. These effects 
serve to considerably enlarge the available parameter space, but in section 4 we show that some 
strong model-independent statements can still be made, and the GMSB models retain a distinct 
character even without taking into account the possibility of discovery modes involving decays 
into the gravitino. 



2. Beyond the minimal model 



In this section we consider a slightly generalized treatment of the minimal model of GMSB. 
The messenger sector consists of a set of chiral superfields ^>j, which transform as a vector- 
like representation of the MSSM gauge group. The supersymmetry breaking mechanism is 
parameterized by a (perhaps not fundamental) chiral superfield S, whose auxiliary component 
F is assumed to acquire a VEV. The messenger fields couple to S according to the superpotential 

W = KS^i'^i (2) 

(Here we have assumed that the messengers obtain their masses only from coupling to a single 
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Figure 1: Contribution to MSSM gaugino masses from messenger field loops 

chiral superfield S; we will comment briefly on the effects of relaxing this assumption below. 
With this assumption, a possible coupling matrix Xij^i^j can always be diagonalized as shown.) 
In the minimal model of GMSB ||2|, $j and consist of chiral superfields transforming as a 
5 + 5 of SU{5) D SU{3)c x SU{2)l x C/(l)y. This choice is sufficient to give masses to all of 
the MSSM scalars and gauginos. 

In the following we will use the same symbol for S and F and for their VEVs. The fermionic 
components of <I>j and <I>i obtain a Dirac mass equal to XiS. Their scalar partners have a (mass)^ 
matrix equal to 

(\^^S? \ 

\\*F* \\S\'') ^"^^ 

with eigenvalues lAj/Spi |Ai-F|. The supersymmetry violation apparent in this spectrum is then 
communicated to the MSSM sector via the ordinary gauge interactions of $i and <I>i. 

The gauginos of the MSSM obtain their masses at one loop from the diagram shown in 
Fig. |l[ The particles in the loop are the messenger fields. Evaluating this graph one finds that 
the MSSM gaugino mass parameters induced are: 

Ma = ^^Y.M^9{x^) (a = 1,2, 3) (4) 

i 

where 



\F/\S^\ (5) 



for each messenger coupling Aj and 



9{x) = ^[(1 + x) log(l + x) + {l-x) log(l - x)] . (6) 

In eq. (|^, na{i) is the Dynkin index for the pair $j,<I>j in a normalization where na = 1 for 
N + N of SU{N). We always use a GUT normalization for ai so that rii = ^Y"^ for each 
messenger pair with weak hypercharge Y = Qem — T3. The variable Xi must lie in the range 
< Xj < 1, with the upper limit coming from the requirement that the lighter scalar messenger 
has positive (mass)^. The minimal 5 + 5 model has = Si'^sC^) = Si^sC^) = 1- Since 
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Figure 2: The functions g{x) and \/ f{x) described in the text. 

the function g[x) obeys g{Q) = 1, in the smah Xi hmit one recovers the result Ma = {aa/4:iT)F/S 
of Q for the minimal model. 

For larger x the expansion 

2 4 fi 

gives good accuracy except near x = I. The function g{x) is graphed in Figure ^, and can be 
seen to increase monotonically with x, reaching a maximum value g{l) = 2 log 2 ~ 1.386. It is 
sometimes convenient to write Ma = j^Aca where 

AGa = ^Y.Mi)9ix^) a = 1,2, 3 (8) 

i 

parameterizes the possible effects of a non-minimal messenger sector and non- negligible Xj. In 
general one finds 

^Na < Aca < 1.386^Na (9) 

depending on Xi, where 

Na = Y.na{i). (10) 
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The scalar masses of the MSSM arise at leading order from 2-loop graphs shown in Figure]^, 
with messenger fields, gauge bosons and gauginos on the internal lines. The calculation of these 
graphs is described in an Appendix, where we obtain the result already given by Dimopoulos, 



Giudice, and Pomarol|18|: 



2|^/^|'E(S) CaJ2 nail) fix.) (11) 



with 

fix) 



r2 



log(l + x) - 2Li2(x/[l + x]) + ^Li2(2x/[1 + x]) 



+ ix^ -x) . (12) 



x^ 

In (|l^), Ca is the quadratic Casimir invariant of the MSSM scalar field in question, in a 
normalization where C3 = 4/3 for color triplets, C2 = 3/4 for SUi2)i doublets, and Ci = 
It is convenient to write fh? = '2J2ai^a/^'^)'^CaA.'g^ with 

Al, = \F/SfJ2naii)fix..). (13) 

i 

In this way the 6 quantities Aca and Asa parameterize the effects of a non-minimal messenger 
sector and non-negligible Xi on the masses of MSSM gauginos and scalars respectively. In the 
limit |F/AjS'^| <^ 1, one recovers the result Aca = Asa = F/S for the minimal model of Q, 
since /(O) = 1. In order to illustrate the relative effects of non- negligible Xi on gaugino and 



sfermion masses, we graph in Figure y the function fix) to compare with gix). When x is 
not very close to 1, one finds excellent precision from the expansion 

/(^) = 1 + ^! _ _ ^l^x^ + . . . . (14) 

■'^ ^ 36 450 11760 ^ ^ 

The function fix) is nearly constant for x not near 1, and falls sharply near x = 1 to a minimum 
value of /(I) = 2 log 2 2 log^ 2 - w 0.702 or V7(T) ~ 0.838. Note that y/7ix) is always 
within one per cent of unity for x < 0.85. Thus as long as \F/XiS'^\ <^ 0.85 for all messenger 
fields, one has simply 

AL = \F/S\'Na (15) 
to a very good approximation. More generally, one finds 

0.838y]V^ \F/S\ iCAsal^VK \F/S\ . (16) 



By combining the bounds on gix) and \/fix) we obtain the result 



^<^<1.65v^ (17) 
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Graph 1 Graph 2 




Graph 7 Graph 8 

Figure 3: Two-loop contributions to MSSM scalar masses involving messenger sector fields. 
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in any model in which ah messenger fields obtain their masses only from a single chiral superfield 
S and its F-term. The efi^ect of non-negligible Xi is always to lower the masses of squarks and 
sleptons compared to the gaugino mass parameters. With some rather mild restrictions, the 
range ( p]7| ) can be significantly tightened. For example, if all Xi < 0.85, one can replace the value 
1.65 by 1.19. With the further restriction that all Xi < 0.5, the same number becomes only 
1.044, so that the scales entering the gaugino and scalar mass formulas differ only at the few per 
cent level. The 1% accuracy level (to which higher-loop corrections are probably comparable 
anyway) for Aca ~ V^^Sa is reached if all Xi < 0.25. 

The masses predicted by equation (^) and (0) are given at the messenger mass scale(s) and 
must be renormalized down to the scale of MSSM sparticles. Decoupling each set of messengers 
^i, at the appropriate XiS, one obtains running DR gaugino mass parameters 

MaiQ) = E MM^^) ■ 



Att S ^ ^ ^ 
\iS>Q 

Below the lightest messenger scale this reduces simply to 

MaiQ) = ^^Ga 



(18) 



(19) 



up to small two-loop corrections [23, 24, 25 1. 



The scalar (mass)^ parameters obtain renormalization group corrections proportional to 
gaugino masses squared, with the result 



rh\Q) = 2j2 Ca 



F 



E 



Qa(Ai5) 



'na{i)f{xi) + 



logQ 



vr 



(20) 



with Ma{Q) given by eq. (0), and aa{Q) by a similar step- function decoupling of messen- 
gers. As long as the couplings Aj do not feature large hierachies, one may safely neglect 
messenger-scale threshold contributions of order 5fh? ~ 2Ca log(Ai/Aj)M^Qa/7r by choosing 
a representative messenger scale Qq ~ \iS. In this approximation one finds 



m 



Q-aiQ) 



Ca 



"a 



where (6i,&2i&3) = (—33/5,-1,3) and 

TaiQ) = [aa{Qo)/aa{Q)f = [1 + (6aaa(Q)/27r) log (Qq/Q)] 



(21) 



(22) 



In the case that all Xi are small and not too different, the running scalar and gaugino masses 
and running gauge couplings can be directly related at any scale by 



(23) 
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while more generally one finds 



= 2Y^CaMl + 1(1 - vM , (24) 

with the ratio ^"sa/^Ga bounded by 0.366/A'a and 1/Na according to eq. (p!7[). These equations 
hold at the one-loop level (with Yukawa couplings and trilinear scalar couplings neglected) 
in a non-decoupling DR scheme, which means that MSSM sparticles and Higgs fields are not 
decoupled at their mass thresholds. In order to make precise predictions about the sparticle 
masses, these parameters must be related to the physical masses of the particles. The necessary 



equations have been given for the gluino and first and second family squarks in |23], and in 



general for all of the MSSM particles in [26| 



So far we have assumed that the messengers all obtain their masses entirely through cou- 
pling to a single chiral superfield S. If this assumption is relaxed, one clearly obtains a much 
more general set of models with a concomitant loss of predictive power. However, the as- 
sumption that only one field S plays a significant role is perhaps sufficiently compelling that 
the alternatives can be considered disfavored. For example, the existence of only one S field 
succesfully addresses the super symmetric CP problem, since all phases in the theory are pro- 
portional to the phase of F/S, and can be rotated away. This need not be so if there is more 
than one field S. The simplest model of this type is perhaps the obvious extension of the 
minimal model of GMSB, i.e. with messenger fields D + D and L + L, and the superpotential 
W = X3S3DD + X2S2LL. The gaugino masses obtained from this model are given by, in the 
small \Fi/\iSf\ limit, 

"3-?^3 „ a2F2 ai /2F3 3F2\ 

^^^ = 4^^' ''' = ^3-2^ = 4^ [53-, ^53-2) ■ 

If the phases of the VEVs are not aligned, this gives rise to an observable CP-violating phase 

arg(MiM|) which could potentially feed into an electric dipole moment for the neutron or 
electron. On the other hand, if squark and slepton masses are very large, such new phases 



could be tolerable, and the interference in (25) could allow Mi to be somewhat suppressed 



relative to M2 and M3 and the slepton masses. However, we will not consider such possibilities 
further here. 



3. Variations in the messenger sector 



One of the outstanding features of the minimal model of GMSB is its predictive power, since 
the values of the soft supersymmetry-breaking MSSM parameters are determined by only a few 
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parameters in the messenger sector. However one can also entertain the possibility of different 
field contents in the messenger sector. The original choice of messenger fields in 5 + 5 of SU (5) 
is motivated by the fact that it is the simplest one which simultaneously provides for plausible 
MSSM masses and maintains the apparent unification of gauge coupling observed at LEP. It is 
well-known that the latter feature is shared by any set of chiral superfields which lie in complete 
SU{5) GUT multiplets. The number of such fields which can be used as messenger fields is 
then limited by the requirement that the MSSM gauge couplings should stay perturbative up 
to the GUT scale Mjj 2 x 10^^ GeV, which amounts to the statement that there can be at 
most four 5 + 5 sets or one 5 + 5 and one 10 + 10. 

While maintaining the apparent unification of gauge coupling is a fine goal, it is not clear 
how much this really should tell us about the messenger sector. First, it is sufficient but 
not necessary to have complete 5 + 5 and 10 + 10 multiplets of SU{5) in order to maintain 
perturbative gauge coupling unification. A counterexample with Ni = N2 = = 3 is a 
messenger sector transforming under SU{3)c x SU{2)l x U{1)y as 

(3, 2, i) + (3, 1, ^) + 2 X (1, 1, 1) + conj. (26) 

which by itself (or with additional gauge singlets) does not happen to form any combination 
of irreducible representations of any simple GUT group. Furthermore, it is not necessary that 
all TeV or messenger scale vectorlike chiral superfields must obtain their masses primarily by 
coupling to the field S. Those that do not can still participate in ensuring perturbative gauge 
coupling unification, but may not act as messenger fields and in particular can have little or no 
effect on the masses of MSSM sparticles. (There is, after all, a precedent already in the MSSM 
of chiral superfields in vectorlike, non-GUT, representations of the MSSM gauge group without 
very large masses, namely the Higgs fields.) Finally, a skeptic might point out that the apparent 
unification of gauge couplings could be partially or wholly accidental, so that it is prudent to 
consider equally all alternatives rather than trust the detailed results of extrapolating coupling 
constant relationships over 13 orders of magnitude in energy. 

Therefore we will consider here the effects of a somewhat less constrained messenger sector. 
We will maintain the constraint that messenger fields should occupy the same representations 
as MSSM chiral superfields. This is motivated by the fact that stable messenger particles with 
exotic SU{3)c x SU{2)l x U{1)y quantum numbers are probably a disaster for cosmology. In 
fact it should be noted that in any case the lightest and the lightest color non-singlet members of 
the messenger sector must be stable insofar as they do not couple to MSSM fields through non- 
gauge interactions. (There is an interesting possibility that a stable neutral messenger might 
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make up the cold dark matter, however |18|| .) Fortunately, small mixings between non-exotic 
messengers and their MSSM counterparts can allow them to decay; the necessary couplings 
may or may not [^] significantly affect the predictions of GMSB. So we consider five possible 
types of messenger fields: 

Q + Q= (3,2,i)+conj.; [/ + U = (3, 1, -^) + conj.; (27) 

D + D=i3,l,^)+coni.; L + L = (1, 2, -^) + conj.; (28) 

£; + E= (1,1,1) + conj. (29) 

with multiplicities denoted {nQ,nu,n£,,nL,nE) respectively. Thus the particle content of the 
messenger sector is specified by a five-tuple of integers, given our assumptions. 

[Actually, as long as we are only using the numbers {nQ,nu,nD,nL,nE) to parameterize 
our ignorance of non-MSSM physics, we can set nu = 0. This is because the gauge interactions 
of any U + [/-type messengers can always be replaced by messengers in the representations 
D + D + E + E, as far as the MSSM sector is concerned, since they have the same index for 
each group. Global features of the theory do depend on nu, of course. One could also consider 
messenger sectors which include single adjoint representations (8, 1,0) or (1,3,0), but we will 
neglect those possibilities here.] 

The number of chiral superfields is limited by requiring that gauge couplings remain per- 
turbative. However, we do not require that the messenger fields by themselves maintain gauge 
coupling unification, for the reasons mentioned above. Instead, we require as our first criterion 
only that the messenger fields should be a subset of some set of fields that maintains pertur- 
bative gauge coupling unification. Assuming that no messenger field mass greatly exceeds 10^ 
TeV, the perturbativity part of the requirement (a^ < 0.2 at Mpianck) amounts to 

Ni = \{nQ+8nu + 2nD + 3nL+6nE)<4: (30) 

N2 = Sng + < 4 (31) 
iVg = 2nQ + nu + nD <4: (32) 

while the full requirement can be written as 

{nQ,nu, riD, nL,nE) < (1, 0, 2, 1, 2) or (1, 1, 1, 1, 1) or (1, 2, 0, 1, 0) or (0, 0, 4, 4, 0). (33) 

It is possible that the requirements (30-3^) can be weakened, but only slightly, by allowing 
the extrapolated gauge couplings to diverge between Mu and Mpianck or by enlarging the 
MSSM gauge group below Mi/. (Additional gauge bosons can contribute negatively to the beta 
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functions for 01^2,3 > but this effect is limited by constraints on proton decay, and by the fact 
that additional chiral superfields which contribute positively to the beta functions must also 
be introduced to break the additional gauge interactions.) The requirements that the gluino 
and the right-handed selectron not be massless at leading order imply N3 > 1 and A'^i > 1/5 
respectively. The possibility A''2 = may not be ruled out yet ||l^, if tan/3 is very small, but 
it should be decisively confronted at LEP2 since it requires a chargino mass smaller than Mw- 
Furthermore, it should be possible to exclude these models with existing Tevatron data if the 
decay A^i 7G is prompt, and perhaps even if it is not. 

There are 66 distinct five-tuples {nQ,nij,nD-,nL,nE) which satisfy these criteria, of which 
53 have 7^ 0. The number of distinct combinations {Ni, N2, N3) arising from these mod- 
els is 40. The ones with 7^ are, in ascending order of A''i: (i,3,2); (|,3, 3); (|,4, 2); 
(1,1,1); (1,3,4); (f,4,3); (1,1,2); (|,3,2); (f,2,l); (|,4,4); (|,1,3); (f,3,3); (2,2,2); (2,4,2); 
(^, 1,1); (^,1,4); (^,3,1); (^,3, 4); (f, 2,3); (f, 4,3); (f, 1,2); (f, 3,2); (f, 2, 4); (f, 4,1); 
(^,4,4); (3,3,3); (f ,4,2); (f,l,l); (f,3,4); (f ,4,3); (f ,1,2); and (4,4,4). The ones with 
N2 = and therefore M2 = at the one loop level are: (|,0, 1); (|,0,2); (|,0,3); (f,0, 1); 
(f,0,4); (f ,0,1); (2,0,2); and (f ,0,2). 

Some indication of the variety which can be obtained is illustrated in Figure which 
shows a scatterplot of Ni/N^ vs. N2/N2, for the 66 models (each shown as an X) which fit the 
criteria (|30|)-(p^. For Xj ^ 1, these quantities are equal to the scale-independent gaugino mass 
ratios {Mi / ai) / {M^ / a^) and {M2 / a2) / {M^ / a^i) respectively. Some of the points on this plot 
are occupied by several models. We have also indicated by circles the presence of 33 models 
which fit the perturbativity requirements (^)-(|3^), but for which (|33|) is not satisfied, so that 
the particle content cannot be embedded into a set which allows perturbative unification of the 
gauge couplings unless additional fields with masses far above the messenger scale are invoked. 

The n X (5 -|- 5) and 10 + 10 models [and the model in eq. (p6|)] all occupy the point 
Ni/N^ = N2/N3 = 1, but there are other models which give quite distinctive and interesting 
predictions. The models on the N2/N2, = axis are the ones with uq = = 0, which 
must have small tan /3 and a chargino lighter than the W boson; we will omit them from the 
discussions to follow. The models close to the N1/N3 = axis have a very large hierarchy 
mg^ <C rriq, and so may be strongly disfavored by naturalness criteria. (We will not attempt 
here a complete analysis of electroweak symmetry breaking requirements.) The most "extreme" 
such model has 

inQ,nu,nD,nL,nE) = {1,0,0,0,0); {Ni, N2, N3) = {^,3,2) (34) 
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Figure 4: A scatterplot of the quantities N1/N3 and for all messenger models satisfying 

the perturbativity constraints (pO|-|32D in the text. Models which do (do not) also satisfy the 
"unification" criteria ( ^3|) are plotted as Xs (circles). In the regime |F/Aj5^| <C 1, the axis 
quantities are equal to the renormalization group scale-independent ratios {Mi / ai) / {M3 / 03) 
and {M2 / 02) / {M3 / a^) respectively. 

with Ml less than gluino and squark masses by perhaps a factor of 50, depending on a^. As can 
be seen already from Fig. ^ there can be quite a wide variety in the mass hierarchies between 
squarks and gluinos and the sleptons and electroweak gauginos. 



4. Constraints on sparticle masses 



In this section we will study some features of the sparticle mass spectrum which follow from 
the 53 models which satisfy the constraints (|30|)-(|33|) and A''^ > as discussed in the previous 
section. We will consider here only the gaugino mass parameters Mi, M2, M3 and the squark 
and slepton masses of the first two families, for which Yukawa interactions can be neglected. We 
will also not concern ourselves with the possible origins or role of the ^ and Bfj, terms or scalar 
trilinear terms. Constraints following from requiring correct electroweak symmetry breaking 
(with viable models for the origins of such terms) will only further tighten the constraints we 
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will derive. In the following we assume F/S < 250 TeV and 0.01 < < 0.99. Taken together 
these imply that the messenger mass scales are bounded by AjS < 2.5 x 10^ TeV. 



It is perhaps easiest to understand the impact of variations in the messenger sector by 
first considering the case that Xi is small for each messenger pair. In that case the quantities 
^i'na{i)g{x,i) and J2i''^a{i)fixi) are each equal to Na, so that the MSSM masses are approxi- 
mately determined by just the parameters F/S, Ni, N2, N3. Using the values listed above one 
can then place some bounds on the ratios of gaugino mass parameters as follows: 

0.067^ < ^ < 3.8^ (35) 
a2 M2 a2 

0.1^ < ^ < 3.4^ (36) 

03 M3 as 

02 M2 .012 . „N 

0.25— < — < 4— (37) 
as Ms as 



Although the gaugino masses run with scale, the veracity of the inequalities (p5[)-(p7|) is renor- 
malization group scale-independent at one loop. [It is not completely accurate, however, to re- 



place ai, a2, as by their measured values at LEP here, since (Pq)-(|37|) hold in the non-decoupling 
DR scheme.] The lower bounds in ( [35| ) and ( |36|) are set by the "extreme" model in (^^. If this 
model is discounted, the values 0.067 and 0.1 are each raised to 0.2. 



The bounds from (|35D-(|37D can be strongly modified by different couplings Aj for messenger 
fields with different gauge quantum numbers. However, some general rules can still be found. 
For example, we find numerically that Mi /Ms is always less than 1 at the gluino mass scale, 
with rough bounds 

iVi Ml A^i 

"■''a5 ^ th ^ "-'as- 

(Note that N1/N3 < 3.4 in these models.) Also, M2 can only exceed Ms at the gluino mass 
scale if N2/N-i > 2, and we always find 

0^21^ < f S 0.6f . (39) 

N3 Ms A^s 

Since N2/N3 has a maximum value of 4 in these models, the overall upper limit is M2/M2, ^ 2.4. 
Similarly, M1/M2 can be as large as about 2.7 at the electroweak scale, when the Xi are chosen 
appropriately and N1/N2 is large. Numerically we find 

"■-t ^ t ^ <-> 

It is interesting to consider the ordering between the mass of the lightest slepton and the 
bino mass parameter Mi, since if is large, this will give an indication whether a slepton or 



15 



a neutralino is the NLSP. Using the approximation of eq. (|2^ ) one finds that 

(41) 



for DR parameters m^„, Mi and ri. Since < l/Ni, one finds that rrigj, > Mi can occur 
only if A'^i < 66ri/(65 — lOri). Now, ri depends on both the messenger scale and the scale at 
which we evaluate the running mass parameters. But a reasonable estimate for the upper bound 
is ri < 1.7 [in the regime of validity of eq. (p^)], from which we learn that mg^ > Mi can only 
occur if A'^i < 2.2. This result still applies in more general situations when eq. (^) must be 
applied. Only 21 models which fit the criteria of the previous section can satisfy this constraint. 
The maximum values of the ratio mg^/Mi in these models are approximately 3.0, 1.7, 1.5, and 
1.35 for (riQ, n[/, n£), n^, n^;) equal to, respectively, (1,0,0,0,0); (1,0,1,0,0); (1,0,0,1,0); and 
(0, 0, 1, 1, 0) (the minimal model). Of course the effect of non-zero Xi can only be to diminish 
the ratio mgj^/Mi, but the electroweak D-term corrections to mg^ can raise this ratio slightly 
if Ml is not too large. There is also a possibility that M2 can be less than both mg^ and Mi, if 
-^1 > -^2- However, even taking into account the effects of non-zero Xi, we find that this only 
occurs for a few models with 

11 11 13 17 19 

(iVi,iV2,A^3) = (y>i>i); (y'1'4); (y-i'i) ^nd (y,i,2). (42) 

These are the models for which a line drawn to the origin on Fig. 4 makes the smallest angle 
with the N2/N3 = axis. 

If Mi,M2 > "mgj^, it is still possible that a neutralino is the NLSP if |/i| is not large. This 
typically means that A'^i has a rather large higgsino content, and A^i jG can be suppressed. 
However, the competing decays A^i hG and A^i ZG may be kinematically forbidden, and 
in any case are subject to very strong kinematic suppressions (1— m|/m^ )^ and (1— m^/m^ )^ 
respectively [^]. Therefore if \/F < 10^ TeV it is still possible to explain the CDF 6677 + 
event with small |//|. This may be particularly plausible in the chargino interpretation ||l2| in 
which the event is due to pp — > CiCi with allowed two-body decays Ci ve and v vNi 
or Ci e^v and Cj, — > eA'^i followed by A'^i — > 7G in each case. Since the production cross 
section for chargino pairs at the Tevatron remains large even for m^^^ ~ 200 GeV, it is sensible 
to suppose that the two-photon event could have been seen even if the decay length of A^i is 
increased by a smaller photino component of A^i. 



The models in eq. ( ^4 ) are also interesting because they minimize the ratio of left-handed 
to right-handed slepton masses. In the regime that all Xi <^ 1, we find that the running 
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mass parameters satisfy m^j^/m^j^ > 1.1 for all of the models which fit our criteria (with 
-^2 ^ 1)- The modification of this ratio due to electroweak D-terms happens to be extremely 
small because of the numerical accident sm"^ 6\y ~ 1/4. However, with appropriately chosen 
Xi, it is possible to obtain mg^ ~ mg^j in the last two models of (p2|). In all other cases, the 
hierarchy mg^ > rn-g^ holds. 

One can similarly analyze the possible ranges for the ratios of squark and gluino masses. 
It is easiest to consider the particular ratio M^^/Mg, since this is least sensitive to electroweak 
effects. Neglecting the quite small effects of U{1)y, one finds for the running mass parameters 

~ o /o 1/2 

^ = ^ [3(A53/Ag3)V3 + (1 - rs)] (43) 

in the approximation of eq. (]2^). This ratio is maximized when N3 = 1 and all « 0, and is 
minimized when A^3 = 4 and all Xj w 1. Thus we find 

^ [1 - 0.67r3]V^ ^ ^ ^ ^ [1 + 2-3]^/^ . (44) 

It is now clear that both the upper and lower limits are saturated when is as large as possible, 
corresponding to a low messenger mass scale. Taking estimated bounds <, 0.72 for = 1 
and rs ^ 0.78 for A''3 = 4, we obtain 

0.65 < ^ 1.48 (running masses at Q = mj^). (45) 



Now, the running masses can be converted into pole masses using the formulas in |23, 26 1, 
yielding a slightly modified estimate for the bounds on the ratio of the physical pole masses: 

0.66 < — ^ < 1.36 (pole masses). (46) 

Mg 



Repeating this type of argument for each value of A''3 separately and taking into account eq. ( P0|) 
one finds 

(0.91, 0.76, 0.70, 0.65) < — <> (1.36, 1.07, 0.95, 0.90) (47) 

Mg 

for the physical mass ratios with A^3 = (1,2,3,4). The upper limits in each case correspond 
to small Xi and small values of XiS. The case of left-handed first and second family squarks is 
slightly different, especially when is relatively large. Numerically, we find 

(0.93, 0.76, 0.70, 0.65) < — ^ < (1.74, 1.23, 1.03, 0.95) (48) 

Mg 

for the physical mass ratios with = (1,2,3,4). 
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The masses of S'?7(2)j;^-singlet squarks are never very different from each other in the models 
of section 3. Taking into account the effects of non-zero Xi, we still find a quite narrow range 

1 < m^^/mj^ < 1.04 . (49) 

This is not surprising since the U{1)y effects are relatively small even for larger A^i. The 
left-handed squarks are always heavier than m-u^ ' "^cZa • Numerically we find 

1 < m^Jm^^ < (1.1, 1.2, 1.3, 1.4) (50) 

for N2 = (1,2,3,4). The squark masses also quite generally exceed slepton masses even for 
models with relatively small N3. Numerically we estimate the bounds 

(1.0, 1.5, 2.0, 2.4) < m^Jm,^ < (4, 6, 8, 10) (51) 

for A^3 = (1,2,3,4). The situation mj^ ~ mg^ only can occur for {Ni, N2, N3) = (^,4,1), the 
highest point in the plot of Fig. 4. 



5. Discussion 



In this paper we have examined some of the possibilities for generalized models of the messenger 
sector of low-energy supersymmetry breaking. Despite the large number of discrete model 
choices and the freedom to vary the Xi = \F/XiS'^\, the parameters of the MSSM are constrained 
in interesting ways. For example: 

• The usual hierarchy mg^ < mg^ < mj^ m^j^ < mn^ mj^ is always preserved, with 
numerical bounds given by (^6|)-(^). 

• The masses of the right handed squarks n^j, dn, Sr and Cr all lie in a narrow band, and 
in a window within about ±35% of the physical gluino mass. The upper limit on mq^/nig is 
(nearly) saturated for the minimal 5 + 5 model with small Xi . The masses of the corresponding 
left-handed squarks can be significantly larger in some models, up to about 1.75m^. 

• The ratios of gaugino mass parameters Mi, M2, M3 can vary quite significantly from the 
predictions of the minimal model, with M2 > M3 and Mi > M2 both possible at the TeV scale. 
However, M1/M3 is always < 1. 

• Only six parameters Aca and Asa [plus the overall messenger scale(s)] enter into the definition 
of the gaugino mass parameters and the first and second family squark and slepton masses. As 
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long as Xi = \F/XiS'^\ is less than about 0.5 (0.25) for all messenger fields, then there are 
only four parameters F/S, Ni, N2, at the 4% (1%) accuracy level, besides a logarithmic 
dependence on the messenger mass scale(s) XiS. 

Let us close by noting a slightly different way to express the constraints on squark and 
slepton masses which follow from the GMSB framework. One can see from the form of eq. ( p4| ) 
that 3 parameters suffice to determine all of the scalar masses for which Yukawa interactions 
can be neglected. This means that for the 7 scalar masses mg^ , m;> , mg^ , m , m^i^ , m^i^ , m 
there must be 4 sum rules which do not depend on the input parameters. Two of these sum 
rules are completely model-independent and should hold in any supersymmetric model (up to 
small radiative corrections [p7|] ) : 

~ = I cos 2/3| ; (52) 
mj^ - ml = M^\ cos 2 (53) 



(We assume tan /? > 1.) The other two sum rules can be written as 



"^L = "^L + 3"^L - si^' cos 2/3| ; (54) 
"^1 = + "^L - \^Ik + si"' ^1^1 cos 2/3| . (55) 



These sum rules are not model-independent. It is interesting to compare with the case of 
models with "super gravity-inspired" boundary conditions featuring a common ttt-q for scalars 
and a common m-]^/2 for gauginos at the GUT or Planck scale. In those models, one finds 



a sum rule which is a particular linear combination of (p4[) and (55): 



2mi^ - - vfi\^ + m\ - mj^ = -y M| sin^ %| cos 2/?| (56) 

This sum rule tests the assumption of a common ttiq. But in GMSB models, one effectively 
has the further bit of information that uiq = (i.e., all contributions to scalar masses are 
proportional to the quadratic Casimir invariants; there is no group-independent piece). This 
leads to the presence of one additional sum rule, which can be taken to be either (M) or (pq). 



It will be an interesting challenge to see to what accuracy these sum rules can be tested 
at future colliders. Perhaps the most interesting possibility is that the sum rules will turn 
out to be violated in some gross way; this would force us to reexamine our assumptions about 
the origin of supersymmetry breaking. As an example, suppose that the messenger sector has 
some feature which causes additional unequal supersymmetry-breaking contributions to the 
diagonal entries in the mass matrix (^). This would lead, through a one- loop graph, to a Fayet- 
Iliopoulos D-term proportional to weak hypercharge manifesting itself in the squark and slepton 
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masses Q. Since such a contribution comes in at one loop earlier in the loop expansion than 
the contributions from the F-term, it is constrained to be quite small in order to avoid negative 
squared masses for some squarks and sleptons. Conversely, even tiny such contributions to 
the matrix (|3|) will be magnified in relative importance, and will therefore quite possibly be 
observable in the sparticle mass spectrum! The impact will be to modify each of the sum rules 
(|5^, (|55|) and (^) by adding contributions — 4L'y/3, 2Z)y/3 and —WDy/S respectively to the 
right-hand sides. 

In general, we find it remarkable that the models discussed here make such a variety 
of testable predictions. In addition to the possibly dramatic collider signatures coming from 
decays of the NLSP into the gravitino, the sparticle spectrum has a rather distinct character. 
Future model building developments will surely tell us even more about what to expect for the 
parameters of the MSSM in the low-energy supersymmetry breaking framework. 
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Appendix 

In this appendix we give details of the calculation of the masses of MSSM scalars which arise 
at leading order from two-loop diagrams involving messenger fields. We employ the component 
field formalism. There are 8 Feynman diagrams which contribute at two loop order, as shown in 
Figure]^. We compute these graphs in the Feynman gauge; then one finds that graph 6 vanishes. 
Each of the other graphs is separately divergent but the sum is finite. It is important to compute 
all gamma-matrix algebra in 4 dimensions before computing the momentum integrals with scalar 
integrands in 4 — 2e dimensions, in order to avoid a spurious non-supersymmetric mismatch 
between the numbers of gaugino and gauge boson degrees of freedom. By straightforward 
methods one finds that the contributions for each pair of messenger fields ^, ^ are given in 
terms of the messenger fermion mass mj = \XiS\ and the two messenger scalar masses = 
\XiS\^ ± \XiF\ by: 

Arh^ = J29aCaSa{^) (sum of graphs) (57) 

a 
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where Ca is the quadratic Casimir invariant [normahzed to {N"^ — 1)/2N for a fundamental 
representation of SU{N)] of the MSSM scalar, and Sa{^) is the Dynkin index of the messenger 
field $ [normalized to 1/2 for a fundamental of SU{N)]. The contributions to the "sum of 
graphs" are given by: 

Graph 1 = -Graph 2 = -^GraphS = 2 (m+) (0,0) +2 (m_) (0,0); (58) 

Graph4 = 4(m+) (0,0)+4(m_) (0,0)- {m+\m+\0) - (m_|m_|0) 

-4m^(m+|m+|0,0) - 4 m^(m_|m_|0, 0); (59) 

Graphs = 8 (m/) (0,0) - 4 (m/|m/|0) + 8 m}(m/|m/|0, 0); (60) 

Graph 6 = 0; Graph 7 = -2 (m+|m_|0); (61) 

Graph8 = 4(m+) (0,0) +4 (m_) (0,0) - 8 (m/) (0, 0) + 4 (m+|m/|0) + 4 (m_|m/|0) 

+4(m^ - mj) (m+|m/|0,0) + 4(m?. - mj) (m_|m/|0,0). (62) 

Here we have used the following notation for euclidean momentum integrals in n = 4 — 2e 
dimensions (omitting in each case a factor /x^*^): 

f d '^q 1 , 

{m,m) = / . . — ^ ^r-TT (64) 

^ ' J (27r)" (g2 + m2)2 ^ ^ 

(milmslmg) = 1 ^ 1 ^^^^^ + ^2)(p + ^2)q^ _ + ^2) (65) 

(mi m2 m3,m3) = / — — - / , , , ^ , 2T772~; 2V7(7 il"; 2^ (66) 

J (27r)"' J (27r)" (g^ + m^)(A;^ + m2)([A; — g]^ + mg)"' 

(In the quantities (0, 0) and (mi|m2|0, 0), it is necessary to keep a finite infrared regulator mass 
nie which will cancel from physical quantities.) 

The terms of the form (m) (0, 0) are easily seen to cancel between the various graphs 1-8, by 
the magic of supersymmetry. The remaining two-loop integrals can be evaluated by standard 
Feynman parameter techniques. First it is convenient to use the identity 

(-1 -|-2e) (mi|m2|0) = mj (m2|0|mi,mi) -|- m| (mi|0|m2,m2) (67) 
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to express everything in terms of dimension 8 integrands. Then one finds, following e.g. the 
methods of l2£ 



(mi|m2|0,0) 



r(l + 2e) 
2(47r)" 



^ + -(1 - 2 log m^) + 1 - vrVe - F2{ml,ml) - 2F3(m?,m2) 



+ [-2 + 2Fi{m{,mi)] logmj + log^ 



+ 0(e) (68) 



and 



(mi|0|m2,m2) 



r(l + 2e) 
2(47r)" 



^ + ^(1 - 2 log ml) + 1 - 7r^/6 - 2 log m2 + log^ m| 



log^ + 2 log mf log m2 — 2Li2 ( 1 — m2 / mf 



2 /™2^ 



where we have introduced yet more notation: 



when a ^ b, and 



Fi{a,b) = {aloga-blogb)/{a-b) 
F2{a,b) = {alog^a-blog'^b)/{a-b) 
Fsia, b) = [aLi2(l - b/a) - bh^il - a/b)]/{a - b) 

Fi{a, a) = 1 + log a 
F2 (a, o) = 2 log a + log^ a 
F^{a,a) = 2. 



+ 0{e) (69) 



(70) 
(71) 
(72) 



(73) 
(74) 
(75) 



The dilogarithm or Spence function is defined by Li2(x) = — ^Q{dt/t)\og{l — xt). Now it is 
straightforward to add all the contributions to the "sum of graphs". In particular, it is easy 
to show that the ultraviolet and infrared divergent terms cancel. The resulting expression can 
be simplified further using standard dilogarithm identities |^^, finally yielding the expression 
given in |18], and in equation (111) of the present paper. 
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